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ABSTRACT 

The Basic definitions for viscoelastic stress intensity factors are set 
forth with associated computations for determining them in laboratory specimens. 
Laboratory tests have been conducted to characterize propellants fracturewise, 
and the most suitable test found was the long strip subjected to a constant 
force loading. The influence of dewetting on the fracture characterization 
appears to be significant, but the matter has not yet been fully resolved. 

An approximate method for performing viscoelastic fracture analyses has 
been presented which incorporates the fore -mentioned fracture characteriza- 
tion. It requires minimum computer usage and should yield good results. 
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VISCOELASTIC ANALYSIS METHOD FOR 



ROCKET MOTORS CONTAINING CRACKS 

INTRODUCTION 

Solid propellant has been found to be susceptible to cracking during 
fabrication or service just like other engineering materials. The cracks 
may or may not be detrimental to the rocket motor depending upon location, 
size, and propagation rate; however a potential danger always exists when 
a crack is present; for instance, the crack introduces additional burning 
which may cause significant overpressure and a subsequent malfunction. It 
is also possible that a crack may create a pathway for burning to the case 
prematurely, producing a burn-through. On the other hand, it may lead to 
other areas of the grain where malfunctioning occurs due to alterations in 
the ballistics or due to loss of structural integrity resulting from the 
high stresses and strains associated with the crack front exceeding the 
design limits. The consequences of such developments must be determined 
for each individual situation; therefore a theoretical basis for predicting 
the initiation and propagation behavior of cracks in propellant is needed. 

The elastic theory of fracture mechanics has been successfully applied 
to rocket motor geometries and loadings by Francis, Lindsey and Parmerter (l). 
Similar work was done on viscoelastic materials by Francis et al (2) using 
an unfilled polyurethane for the grain material. Two-dimensional grain 
cross-sectional slices containing cracks were analyzed and compared with 
experiment. Predictions of viscoelastic crack behavior were made by using 
a sequence of elastic crack propagation problems. The method of fracture 
characterization followed that used by Knauss(3) on the same polyurethane 
material. Progress was made and results of the study were used to 
provide guidelines for some of the concepts developed in this work. 
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In this report the fundamental expressions for stress and displacement 
near a crack tip are developed for linear, unfilled, viscoelastic materials. 
These expressions are used to reduce crack propagation data obtained on strips 
of propellant that are considered to be infinite in extent. Both constant 
stress and constant displacement boundary conditions are employed in the tests 
for purposes of comparison. Experimental fracture characterization in terms 
of stress intensity factors versus crack velocity in propellant is the end 
result. Finally a method is presented for the viscoelastic analysis of rocket 
motors containing initial cracks. 
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DEFINITIONS AND ASSUMPTIONS 



Stress Intensity Factor 

Elastic solutions for the stress field in the vicinity of a crack by 

(4) (5) (6) 

Williams , Westergaard ' and Irwin ' have been written in terms of 
undetermined coefficients that became explicit when the particular body 
geometry and loading are specified. In terms of the mixed coordinates 
(see Fig. 1) 




Figure 1 Stationary crack coordinate system. 



normally used to express the stress field, 



K 1 

a.,(r,e) = -±- 

J J 2r 



f. .(0) + 
ij 



K o 

— = s^O) 

J 2r J 



( 1 ) 



where and are the opening mode and shearing mode stress intensity 
factors (S.I.F.) respectively. 

To more precisely define these quantities, use can be made of the 
opening and shearing stresses in expanded form. 
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These two stress components possess a certain assymetry, and the S.I.F. can 
he conveniently defined in terms of them. 

K = lim V§ra (r,0) (3a) 

r 0 y 

K 2 * r 1 " m o^ ? V <r ’ 0) (3b) 

where the limit is used to indicate that this relationship only applies in 
a small neighborhood of the crack -tip. 

For 9 =0, r = x f and y = 0. Making this substitution, the definition 
can be written also in terms of cartesian coordinates. 



K 



= lim J2x' o (x',0) 
x - 0 * 



K = lim ^/2x ' t (x',0) 

2 x - c> 



(4a) 

(4b) 



In viscoelastic fracture no standard definitions of stress intensity 
factor, fracture toughness, coordinate orientation, etc. have emerged; 
however it appears straight forward to extend the above definitions in the 
obvious way. Since the crack is normally propagating under load in visco- 
elastic materials, it is convenient to shift the base of reference to the system 




Figure 2 Moving crack coordinate system. 



Translating the previous expressions for S.I.F. by an amount c(t) and 



allowing for the possibility of the stresses varying with time, 

K v (t)= lim {[2(x - c(t))] 2 CT (x,0,t)} (5a) 

1 x - c (t ) y 

K 2 V (t) = lim {[2(x - c(t))] 2 x (x,0,t)} (5b) 

x - c ( t ) ^ 



Apart from the square -root- of -two factor, this definition is the same 
as used by Graham^ \ Knauss^^ also used this form, but Meuller^^ and 
Knauss and Meuller^^used a nondime ns ional form of the above for an 
infinite strip of width 2b. 



I = lim 
x -* c(t) 



r p CT y (x?0,t) ' 



(6) 



Since all of these forms are closely related, and since the definitions 
of equation (5) possess the appropriate limit check for elastic materials, 
they will be used throughout this report. 

Fracture Postulate 

Fracture studies in recent years have led to the conclusion that flaws 
in materials will grow when under load. This continues until the flaw 
becomes of critical size, at which time the specimen undergoes unstable 
crack propagation and fracture . This has been found to be true for metals , 
especially in fatigue, and for nonmetals alike. 

(O \ 

Knauss' has published theoretical developments of these ideas for 
viscoelastic materials and has shown that the important parameter of 
interest for rate sensitive materials is crack velocity. In the elastic 
case, the threshold of rapid crack propagation is often the focus of 
interest; however in viscoelastic materials the crack is found to propagate 
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for almost any load level. The velocities are infinitesimally small for 
low loads and increase as the body is more highly stressed. The question 
to be answered then becomes , not at what load level does the crack propa- 
gate, but at what velocity does the crack propagate for any given load. 

As defined, the viscoelastic stress intensity factor contains information 
relative to the specimen geometry, material and boundary conditions; 
therefore we desire the relationship between stress intensity factor and 
crack velocity. 

This brings us to the fracture postulate: 



In words this states that for every stress intensity factor there is a 
corresponding unique crack velocity. 



have duplicated these tests and done additional work on other types of 
loads.) For all temperature ranges tested, the material was elastic in the 
sense that it did not relax; however the crack behavior was viscoelastic 
and exhibited a velocity that varied with loading. Knauss postulates that 
in this type of material, when the strain at the crack tip is the same in 
two different bodies, the same crack velocity will be produced. Of course 
in a body that is behaving elastically so far as stress and strain are 
concerned, the same strain level is equivalent to the same stress intensity 
factor. Thus, in the postulate of fracture, we are following the previous 
ideas of Knauss. Furthermore, it states that if the stress intensity factor 
history can be calculated, the crack velocity equation can be determined 



If 




then 




Knauss has performed experiments on Solithane 113, a polyurethane 

( 2 ) 

elastomer, which have substantiated this to a degree (Francis et al v ' 



(3) 
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from the experimental characterization of the material. Once this is 
known, the crack history for a rocket motor can be obtained. This will 
be demonstrated later. 

EXACT SOLUTIONS FOR S.I.F. 



The classical correspondence principle normally applied in the solu- 
tion of boundary value problems in the linear, quasi-static theory of 
viscoelasticity, does not apply for crack propagation problems. In this 
instance new boundary is being created as the crack propagates , and the 
conditions on the boundary are changing in type; i.e. at the crack tips 
they are changing from displacement to stress boundary conditions and vice 
versa. In this kind of problem there will be points of the boundary where 
only partial, histories of the boundary conditions are known a priori; 
consequently the integral Laplace Transform cannot be computed, and the 
correspondence principle is not applicable. 

Graham^*"*" ^ has extended the correspondence principle to encompass 
crack problems . Two examples that are of particular interest in fracture 
characterization of propellant are the long strip with either displacement 
or stress boundary conditions . 

Long Strip with Stress Boundary Conditions 

The sheet and loading are symmetric about the x axis , as shown in 
Fig. 3, where the lower half of the stip is removed and replaced by the 
stresses which act along the sheet's horizontal centerline. 
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Let the entire boundary of the sheet be defined as (3 with p^t) 
representing the uncracked portion of the x-axis and P-^t) representing 
the remaining boundary. Then p = P 1 (t) U P 2 (t). The boundary conditions 
are : 



a s (X,t) = 


0 






on 


p 


1 


f° 


on 


crack edge j 






CT n (X,t) =< 


0 


on 


x = 0,L \ 


on 


p x (t) 




^ a (x,t) 


on 


y = b J 






u n (X,t) = 


0 






on 


p 2 (t) 



(7) 



the subscripts s and n respectively refer to normal and shear components. 
The boundary conditions of Equation (7) meet the criteria of Graham's 
extended correspondence principle. The extended correspondence principle 
permits the construction of the viscoelastic solution for a moving crack 
from an elastic solution for a moving crack. 



a . V (t) = L' 1 [ a . 1 ? 6 * ( s ) ] E -> E* 

o - u* 



(8) 



where m designates a moving crack. The moving crack solution may be 
generated from a similar stationary elastic crack solution. For an 
incompressible elastic strip, the stresses in the neighborhood of the crack 
tip are given by, 



e = f (Q) 



(9) 



J 2rn 

where a = boundary stress in uniform stress region 
2b = strip width 

The crack is now imagined to propagate slowly. For most elastomers 
and all propellants, the velocity of propagation is quite low, which makes 
the inertial terms in the equations of motion small and negligible. The 
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field equations become the same as those of the static crack problem; 
therefore the stress field is the same as that given in Equation (9) , 
except that the coordinate system to which the above quantities are 
referenced is now moving since it is attached to the crack tip. The 
dynamic, elastic solution is 



g? (t) = jzMZH f re(t)i (10) 

J V2nr(t) 1J 

It is a simple matter to reference these quantities to an inertial 
coordinate system. As stated previously the loading and geometry is taken 
such that only the opening mode of propagation is present; therefore from 
Fig. 2 

r(t) ' cos Q (t) + c(t) =x (ll) 

Equation (10) for the dynamic elastic problem becomes. 



me (t) , F [e(t) ] 

J / /2tt(x - c(t) 



( 12 ) 



To use the correspondence principle, the transform of Equation (12) is 
needed. 



me* _ r q(t) /b~ 



me* r 
a . . = 
ij L 



J 2n[x - c(t )] 






(13) 



By Equation (8), the corresponding viscoelastic expression is 

ctT* = f F..,re(t)if (14) 

1J V2n[x -c(t)] 1J J E -* E* 

V “* V* 

Since E and v cLo not appear, the inversion is trivial, and the viscoelastic 
stresses are identical to the elastic ones. 
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(15) 



CT 




CT ( t ) F. .[© (t)3 

J 2tt(x - c(t)) 1J 



From this result and by application of the definition of stress intensity 
factor given in Equation (5) 3 



k: 



(t) 




(16) 



This viscoelastic solution is accurate to the degree that the elastic 

solution from which it was constructed is accurate. In this case it is known 

(l2 ) 

to he good via Knauss series solution. 

Long Strip with Displacement Boundary Conditions 

Many elastic crack propagation studies have been conducted on strips 
with displacement boundary conditions, but in order to use it for viscoelastic 
fracture, the viscoelastic solution is needed. This is not directly 
obtainable from Graham, but it is possible to gene ate it by following some of 
Graham’s methods. 

Consider the viscoelastic boundary value problem of Fig. 4, where A, B, 
and C do not involve material constants, and the crack is propagating. 




Fig. 4 Viscoelastic problem. 
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Following Graham, we consider a dynamic elastic problem with a moving crack 



and the same boundary conditions. This is shown in Fig. 5. 




Fig. 5 Elastic problem. 



0 0 

The solution to this problem, a. .(x,t) and u.(x,t) is known, and with 

i J l 

it, all field quantities can be computed on the boundary. This solution 
allows us to define the equivalent elastic problem shown in Fig. 6. 



u e = 0 
n 



e 

x 

s 




u = 
n 

e 

u = 



0 

C(x,t) 



Fig. 6 Equivalent elastic problem. 

The transform of this problem, which can now be computed because the boundary 
conditions of the moving crack are known for all time, is shown in Fig. 7. 
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•* 

u e = 0 
n 

* 

0 * 

u = C 
s 



n 



Fig. 7 Transform of equivalent elastic problem. 



Now we define an equivalent viscoelastic problem whose boundary 
conditions are defined .in a manner to "correspond" with Fig. 7. This is 
shown in Fig. 8. 




u 



= 0 



< = C(x,t) 



CT 



V 

n 



= 0 



Fig. 8 Equivalent viscoelastic problem. 



In Fig. 8, K(t) is defined as 



K(t) - L^j 



(17) 
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The transform of this problem can also be computed by virtue of the 
definition of boundary values for all time. The results are shown in Fig. 9 . 




v* 



u, 



n 



v* 



u 



s 



0 

C* 



Fig. 9 Transform of equivalent viscoelastic problem. 



Comparing the governing equations and boundary conditions for the problems 
shown in Fig. 7 and Fig. 9? it is seen that a correspondence exists. 



v*, ^ r e *"l 

= L'UjJe 



E* 



(18) 



This will also be a correspondence principle for the original problem 
of Fig. 4 if it can be shown that on the extension in front of the crack 
the boundary condition of 



°n = f dT («) 
o- 

produces u^ = 0 over the same surface. It goes without saying that the zero 
shear stress condition must exist due to the assumed symmetry of loading and 
geometry. 
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Using the extended correspondence principle, the displacements on 



the boundary are 



C<*.P> = [ u n* (x ’ p) ]E~E* 



( 20 ) 



e e* 

by symmetry u =0, and u =0. If there are no material constants in 

e* / \ 

, then from (20), 




( 21 ) 



and on the boundary 




( 22 ) 



Thus, the problem shown in Fig. 8 is the same as the original one in Fig. k 
for which a solution is being sought, and a correspondence principle has 
been established for this case. 

Applying Equation (l8) to the elastic stress solution of the long 
strip with displacement boundary conditions , we obtain 




(23) 



Inverting the transform 




(24) 



where 




If we look at a point on the line of propagation 



e(t) = rr h(t - t 1 ) 



(25) 
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where is the time of arrival of the crack at the point of interest. 

Using Equation (25), it can be shown that f in Equation (24) reduces to 

yy 

f yy[ 0(t) ] = h(t) ' h(t ' V (26) 

For a step displacement input, v Q (t) = V Q h(t), the expression for opening 
mode stress simplifies to. 




= V c 4j 



o- 



J REL 




J2 (x - c(t) ) -I 



(27) 



Carrying out the differentiation, and simplifying 



v , v ,k_ f E KEL^ t ^ p t E REL^ ~ 

’ " ° I,, 1 [2(x . c(t)]3/2 



dT 



] 



( 28 ) 



With this expression for hoop stress, the stress intensity factor can 
be obtained via the definition of Equation (5). 



V/ s lim ^ r— t Ik r E REL (t) . P t n “reL^ ' T;c '' T '' , -] 

K ± (t) ■ X „ “ (t) SzUTTTm v o ♦ J i [g(x . c(t))]3 / ^] (29) 



,t, E REL (t ' T) ° (T> 



It goes without saying that this equation is difficult to use in data 
reduction for crack propagation tests. This has led to a search for an 
approximate solution possessing a simpler, more usable form. 

APPROXIMATE VISCOELASTIC STRESS DISTRIBUTION 



The expressions for the elastic stress distribution in the neighborhood 
of a crack tip can be written as 



K® (E,a,c) 



is 



ijs 



V2r 



f. .(e) 

ij 



(30) 
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where superscript e denotes elastic and subscript s denotes stationary crack. 
The functional dependence of K upon E ,a, and c, indicates that the stress 
intensity factor is in general a function of material, loading and crack 
length. Because the crack is stationary, the viscoelastic correspondence 
principle is applicable. 



a 



v* 

ids 



K.®*(E*, a *,c) 



J 2 



f. .(0) 
ij 



(31) 



Inverting the transform, 



a 



v 

ijs 



-ZZ. (E(t),a(t),c) f..(0) 

J2r 1J 



(32) 



where as before, to fit the definition, 




(33) 



For slowly moving cracks , the inertia terms in the equilibrium equations 
can be neglected, and the field equations for the moving crack become 
identical to those for the stationary crack If it were not for the changing 
boundary conditions at the crack tip, the t > solutions would be identical; 
however the S.I.F. is altered as the crack propagates because the geometry of 
the body is changing and 

■l. K (3-) 

In the case of the long strip laboratory sample, the specimen approxi- 
mates on infinite strip containing a semi-infinite crack. As the crack 
propagates, the geometry remains unchanged to the degree that the specimen 
models the infinite sheet; in fact, any geometry for which the crack propa- 
gation makes minor perturbations on the nature of the body, such as infinite 

domains , will be well represented by 
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(35) 



K^(t) = [K 7 ] v 

l v 7 L 1S J C - c(t) 

To illustrate this , consider an infinite sheet with a central crack of 
length 2c. 

K is = cr ^ T 06) 

By the correspondence principle and the definition of Equation (33) , 

K is = L_1 1°*^ = (37) 

Using the approximation of Equation (35) > the viscoelastic S.I.F. for a 
moving crack is found to be 

= CT(t)VcTty (38) 



In this case, this is exactly the expression that results from the more 
rigorous extended correspondence principle. 

Infinite Strip With Displacement Boundary Conditions 

Beginning with the elastic solution for a stationary crack in a long 
strip , 



a 



e 

ijs 



■4jT 



v o (t) 



3bn \/STrr^7 J 



f. .(e) 



(39) 



By the classical correspondence principle, 





f 



V * 

o 



V2(x - c) 




(40) 



Inverting the transform, still holding the crack stationary, 



v = J£ *V 9) .- r* E 

J o- 



dv 

REL^ " 17" 



dT 



(41) 
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For a step displacement loading, 



v = V h(t) 
o o v 7 



V 

a. . = 



Ik F i;< e) 



V E Drr (t) 



1Js 3bn ^TT7T 0 EEL 



(42) 



Converting to a moving crack through the approximate relation of Equation (35), 
9(t) is defined as before, 



ct V = Erel ^^ 

CTid 3bn o y 2 ('x - c(t) 



[h(t) - h(t - t 1 )] 



( 43 ) 



The stress intensity factor is readily recognizable from this form as 



K I (t) - ^ V o E EH, (t) 



( 44 ) 



In order to measure this quantity on each specimen directly which would 

account for any specimen to specimen variation in E n _ T (t), we seek to express 

KEL 

the S.I.F. in terms of stress. For a long strip without a crack, the princi- 
pal elastic stress is 



V 

e _ E o 

CT y 2 N b 



(1-u ) 



( 45 ) 



Applying the correspondence principle to an incompressible strip, 



a = 4i- E* V * 

y 3b o 



( 46 ) 



Upon inversion, 



dv 



CT y^ 3b J_ e rel dT dT 



O- 



( 47 ) 
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For a step strain history, v q = V Q h(t) 



<V (t > - sr v o 



(48) 



From Equation ( 44 ) 



* 1 ^ - <£ <*> 



( 49 ) 



By using the stress expression for an uncracked strip , a V (t) represents 

y 

the stress in the undisturbed field far out in front of the crack. In order 
to reduce laboratory data, it is necessary to relate the total load on the 
specimen, or the displacement of the specimen grips, to the far field stress. 
This has been accomplished through a finite element program. 

v- . 

Because the stress, a (t), is constantly relaxing during the test, the 

y 

stress intensity factor is likewise decreasing. This results in the crack 
velocity steadily diminishing and approaching zero quite rapidly. The 
constantly varying crack velocity and, more particularly, the fact that it is 
approaching zero makes this test somewhat unattractive; however it is not 
impossible to use. 

Infinite Strip with Stress Boundary Conditions 

The same method when applied to the infinite strip with stress boundary 
conditions leads to considerably different results. Beginning with the 
elastic solution, (ll) 






( 50 ) 



The viscoelastic stress intensity factor for a stationary crack can be 
obtained from Equation ( 33 ) 




(51) 
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Since the stress intensity factor is independent of crack length, it is the 
same for a slowly moving crack as for a stationary one. 



This agrees exactly with the rigorous solution obtained in Equation (l6). 

Thus if stress is held constant, the stress intensity factor remains constant 
and so does the velocity. However, this test is very difficult to run in 
actual practice. Because the specimen is not actually infinite, the net 
cross-section stress is difficult to hold constant as the crack runs. 

( 13 ) 

Hertzler has painstakingly done this in order to obtain data for this 
particular loading condition; however it is much more convenient to work with 
constant load. 

Infinite Strip With Constant Load Boundary Conditions 

Probably the most convenient approach to viscoelastic fracture character- 
ization from the experimental point of view is to apply a constant load. A 
strip specimen is bonded between rigid grips as in the case of displacement 
boundary conditions. (See Figure 10.) 




(52) 



f =. Constant 



<r 




<r 



F 



Fig 10. Constant Force Sheet 



20 



The S.I.F. for this specimen in terms of the far-field stress is given in 
equation (49)* 

K^(t) = o(t) J ( 53 ) 

The relationship between o(t) and the applied force has been obtained from 
a finite element solution for a 1 " x 6 " specimen which has been used in 
this research program. For crack lengths up to three inches, the difference 
between the far-field stress in the undisturbed region and the average stress 
based on net cross-section was less than 1 $. Considering this error to be 
negligible , 



=j$r 




( 54 ) 



where F(t) is the applied force and A(t) is the net cross-sectional area 
computed from the specimen length minus the crack length. 

The variation of A with crack length will be very slight when the crack 
is small compared to the sheet length, which should be the domain of most, 
if not ail, of the test. By varying the level of applied force, the velocity 
of crack propagation in the sheet is varied. For a given force, the S.I.F. 
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is virtually constant as mentioned before; consequently the velocity of 
crack propagation will also be constant. A plot of vs. c(t) completes 
the characterization . 

Infinite Strip-Constant Displacement Rate 

A good test case to which the fracture characterization might be applied 
is the infinite strip subjected to a constant displacement rate loading. 

v 

Having obtained the fracture characterization of vs. 6 via the constant 
load test, it is possible to predict the crack behavior for the same strip 
under constant displacement rate and check the results experimentally. 

From the elastic solution. 




Using the correspondence relation of Equation (33)? 



(55) 



K 



.V* 



~ E*v * 
3brr o 



Inverting the transform, 



Kj(t) = 



/T P t w , + , . 

'3hJ E REL ^ ” T d T T 
o- 



For constant displacement rate, R, 



(56) 



(57) 




(58) 



The integral can be evaluated by means of a Prony Series representation of 
the relaxation modulus. 



E REL^^ 



E.e 

1 



-c^t 



+ E 



(59) 
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Substituting into the above 




-a-tvn 
1-e 1 I 



)j + V} 



(60) 



From this expression and the fracture characterization, a crack history can 
be predicted. This solution is applicable for crack lengths that approximate 
the infinite strip condition. Once the crack becomes great enough, by 
comparison to the length of the strip, to cause the far -field stress to deviate 
from the net section stress when it runs, the above expression will begin to be 
in error. The amount of the error can be determined numerically. 



The first two methods of fracture characterization previously discussed 



be carried out in the near future. 

Constant Displacement Case 

Tests were conducted on cracked biaxial strip specimens with displacement 

boundary conditions. The strips were six inches in length, one inch in width 

(2b) and l/lO inch thick. Redwood tabs bonded to the horizontal edges allowed 

for gripping. An initial crack of 0.75 inch was cut, as suggested by 
(9) 

Mueller ' , and the length of the propagating crack was measured as a function 
of time with an optical comparator. 

The testing was performed with the apparatus, which consisted of an 
Instron Universal Testing Machine equipped with special jaws to clamp the 
specimen and apply displacement v q while maintaining the upper and lower 
boundaries parallel. The displacement was applied at rates ranging from one 
inch/minute for low temperatures to five inches/minute at the high tempera- 
tures. A conditioning chamber allowed the tests to be conducted at the same 
specific temperatures as in the previous case. After waiting ten times the 
loading time the ramp tests simulated step displacement loading, v Q h(t), for 
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which the stress intensity factor reduces to Equation (44). As time 
proceeds after loading, the stress intensity factor is monotonically 
decreasing as a result of relaxation. This produces a corresponding decrease 
in crack velocity, which was observed experimentally. 

The stress intensity factor for each test was calculated by using 
Equation (49) and plotted on a graph of stress intensity factor versus crack 
velocity for each specific testing temperature (see Figs. 11 through 16). A 
conventional a^ shift was attempted with the data, and a smooth master curve 
was obtained graphically by shifting the curves for various temperatures along 
the abscissa, relative to the curve at T q to produce a continuous function. 

The resulting characterization curve is shown in Fig. 17 where T q = 70°F, 
and a^ = 1. It turned out that the a^ curve for the fracture characteriza- 
tion corresponded exactly with the a^, curve for modulus. 

Constant Stress 

The specimen configuration used to achieve the stress boundary conditions 
is shown in Fig. l8. The specimen was cut from samples designed for use in 
biaxial tests with a gage length of 1.0 inch and thickness of 0.1 inch. The 
wooden tabs were segmented with saw cuts, which stopped in the adhesive bond 
just short of the propellant . (To penetrate through the bond into the 
propellant would have created a stress riser and a point of probable crack 
initiation.) The width, w, ideally would be very small, allowing shear-free 
edges; however, practical considerations forced a width of w = 0.25 inch. 

Initial experiments, with the right-hand edge of the propellant flush 
with the wooden tab edge, showed that soon after loading a fast running crack 
would initiate at a local defect along the right edge of even a carefully cut 
specimen. Adding a l/8 inch unloaded length reduced the stress along this 
edge and eliminated the problem. To compensate for this departure from 
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FIGURE 12 PROPELLANT FRACTURE CHARACTERIZATION CURVE FOR 30 
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FIGURE 15 PROPELLANT FRACTURE CHARACTERIZATION CURVE FOR 50 °F 
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FIGURE 14 PROPELLANT FRACTURE CHARACTERIZATION CURVE FOR 70 °F 
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FIGURE 15 PROPELLANT FRACTURE CHARACTERIZATION CURVE FOR 75 



Stress case 
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FIGURE 16 PROPELLANT FRACTURE CHARACTERIZATION CURVE FOR 85 
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FIGURE 17 PROPELLANT FRACTURE CHARACTERIZATION MASTER CURVE 
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Fig. 18 Specimen configuration for stress boundary conditions. 
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uniform stress at the right hand edge and yet still meet the restriction 

leading to Equation (52), c(t) was limited to 2.0 inches. 

max 

The thickness of the saw cuts in the wooden tabs was large enough to 

allow unrestrained e , or d/w ^ e . The upper segments were individually 

x x max 

suspended and the lower segments were individually loaded by weights W. An 
initial crack of one inch was cut with a razor blade. Horizontal extension 
of the crack and constant propagation velocity were taken as evidence of 
uniform loading. Each specimen was tested by applying all loads at t = 0, 
and the length of the propagating crack, c(t), was measured with an optical 
comparator. As the crack passed a segment, the weight on that segment was 
removed in proportion to the width of the segment traversed by the crack. 

The tests were conducted in an environmental chamber with the relative h umi dity 
less than 30$> at six temperatures ranging from 15°F to 85°F. The crack 
velocity was slower in this geometry for a given S.I.F. than it was under 
constant displacement as is shown in Figs. 11 through l6. The shifted master 
curve for constant stress, which is shown in Fig. 17, was obtained by using 
the same a^, shift factor as was determined in the first set of tests for 
constant displacement. 

Differences in the fracture behavior of the load histories is thought to 
be due to dewetting. The theory developed considers only a purely 
viscoelastic material. However, the solid propellant was a filled material, 

(14) 

which when loaded, exhibited the additional mechanism of dewetting. Wood 
has shown that propellant dewetting is not a rate process, rather it occurs 
upon initial application of stress. This behavior has a significant bearing 
on the experimental results obtained in this study as will new be shown. 

In both cases tested, the measured load parameter was stress, thus 
automatically accounting for the modulus relaxation in the displacement case. 
However, the two stress histories were quite different as shown in Fig. 19. 
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In the stress case, stress was constant while in the displacement case there 
was an initial extremely high stress which continually relaxed. 



or <r 




Constant Stress Constant Displacement 



Fig. 19 Stress histories. 

Thus a sample tested in the stress case exhibited min i m um dewetting 
because of its history of constant stress. However, in the displacement 
case, for a corresponding stress at t 2 lOt^ ( w ^ ere ^ me re< l u i ret ^ 

to apply the load) much more dewetting had occurred due to the initial high 
stress. Consequently the crack propagated through a material which had 
been significantly altered from that in the stress case. The propagating 
crack had fewer bonds to break per unit area, compared with the stress case, 
and accordingly its velocity was as much as two orders of magnitude higher. 
This produced significantly different characterization curves for the two 
cases as shown in Fig. 17. 



Since dewetting produces microscopic voids in the material, stiffness 
is reduced as dewetting increases. Using this fact, the increased dewetting 
of the displacement case was verified in preliminary tests at 15°F. In 
these tests, for identical stress and time, the modulus in the displacement 
case was one-half that of the stress case. 

In an attempt to reduce the initial high stress peak in the displacement 
case, thus reducing the dewetting, additional preliminary tests at 15°F were 
conducted. The maximum stress was limited to that which produced stress 
intensity factors comparable with those of the stress case for the same 
temperature. Equivalent moduli showed that the amount of dewetting was in 
fact equivalent. However, as soon as the displacement rate was stopped, the 
stress relaxed due to Thus even if the crack propagated at the same 

velocity as in the stress case (approximately 0.1 inch per hour) relaxation 
of the driving stress decelerated the crack to zero velocity over distances 
too small to be accurately measured. Thus, no additional characterization 
data were obtainable with this loading method. 

The accuracy of the measured velocity in the stress case was considered 
superior to that of the displacement case, especially for low velocities. 

For in the displacement case, c(t) had to be measured over relatively short 
periods of time to obtain an accurate crack propagation function. (Velocity 
was obtained by graphically measuring the slope of this function.) As c(t) 
became small, Ac(t) for the time interval mentioned above also became small. 
Thus, because the experimental crack propagation was irregular due to 
inhomogeneities in the propellant, the actual velocity was obscured as those 
irregularities approached the magnitude of Ac(t). However, in the stress 
case, c(t) was a constant and Ac(t) could be measured over large time intervals. 

The scatter in the characterization curves is due to the inaccuracies 
mentioned above for the displacement case and also to the inhomogeneous nature 
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of a highly filled viscoelastic material such as propellant. By compari- 
son, similar characterization of solithane, an unfilled and relatively 

( 9 ) 

more homogeneous material, by Francis et. al. ' showed significantly less 
* 

scatter. 

Referring to the characterization curves , in both cases there appears 
to be a critical below which a crack will not propagate. Again for 
both cases, especially apparent in the stress case, there appears to be a 
maximum velocity of propagation; however these maximum velocities are 
somewhat misleading. Preliminary tests showed that at very high loads the 
failure was general and catastropic, initiated at voids and other stress 
risers throughout the specimen as well as at the crack tip. Clean crack 
propagation was possible at high loads for low temperature, but high loads 
at room temperature precipitated fracture throughout the specimen 
resembling disintegration. 



*Solithane at room temperature has a relaxation time of less than one 
minute. Thus, though tested for the displacement case, essentially constant 
velocity crack propagation in these tests allowed the additional advantage 
of increased accuracy of the measured crack propagation function for reasons 
mentioned above. 
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ROCKET MOTOR ANALYSIS 



The purpose of all of the preceding material relative to making a 
fracture characterization is to obtain the base information required to 
analyze a rocket motor. Rocket motor geometries are generally complex 
and finite element methods are required to generate solutions and stress 
intensity factors even for elastic fracture predictions. A method will 
now be presented which uses the elastic methods to generate the desired 
viscoelastic fracture results. 

Viscoelastic fracture analyses of rocket motor geometries compounds 
the two complexities of time dependent stress analysis and fracture 
mechanics. The expense of the undertaking, as well as the difficulty of 
doing it rigorously, suggests that approximate methods of solution be 
pursued and developed. 

Appendix I has been written to give motivation to an approximate 
method of stress analysis commonly referred to as quasi-viscoelastic. In 
this development, a circular cylinder is used as the motor geometry, and 
it is assumed that the method holds for all geometries. With a means 
available for obtaining the stress field, another approximate analysis 
procedure is now presented to determine fracture behavior. 

For most, if not all, two-dimensional geometries with through-the- 
thickness cracks, the stress intensity factor has a common form. 

K 1 = \a /c~ ( 61 ) 

where \ is a modifying factor containing information relative to geometry, 
and in a sense, reflects the departure of the geometry in question from 
the infinite sheet with a central crack where X = 1. X includes 
influences such as boundaries, stiffeners, dimensionality of the stress 
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field, etc. a is the applied uniaxial stress, and when cast in this form 
it is seen that represents the stress that would exist in the body at 
the location of the crack if there were no crack present. In other words, 
o is the stress acting across the line of the crack in a identical 
geometry but without a crack. For the rocket motor shown in Fig. 20, 




Fig. 20 

the crack extends 

crack is . The 
0 

K x = \o 0 J~T (62) 



along the r direction, and the stress acting across the 
stress intensity factor for the motor can be written as 



where the crack length is c, and the influence of the bore, case, trans- 
verse stress and material properties are included in X . For a given 
configuration; the functional dependence of the quantities in Eq. (62) can 
be represented in terms of these parameters which will vary with time in a 
viscoelastic fracture analysis. 

K le = X(E,c) ct 0 (E,p) JT (63) 
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The subscript, e, designates elastic stress intensity factor, while E, c 
and p represent the modulus, crack length and pressure loading respectively. 
These are the quantities that will vary with time in a viscoelastic fracture 
analysis . 

It is proposed that \(E,c) be found by numerical methods using finite 

( 2 ) 

element codes. Such a study was conducted by Francis' , et. al. for 
elastic materials and they experimentally confirmed predictions made in 
this manner. For the geometry specified, the elastic stress intensity 
factor was computed from energy differences for runs with incrementally 
different crack lengths . 



2 = E AU 

le Trt A c 



(64) 



( 15) 

(The procedure is defined in a note by Deverall and Lindsey' .) Results 

are shown in Fig. 21 for a high modulus material. Similar curves could be 

constructed for the range of modulus values exhibited by propellant, and X 

follows by dividing K, by a reevaluated at r = a + c. 

j-0 y 




Fig. 21 

Stress Intensity Factor vs. Crack Length 
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Two examples for low modulus materials were available from the 
literature in reference 2 for pressure loads. They are reproduced in 
Figs. 22a and 22b, and from them X was computed to show the dependence of 
this factor on modulus. Of course this range of modulus is much greater 
than that found in rocket motors. However once these results are obtained 
X is easily computed. 

A plot of the results in Fig. 23 shows that X is not strongly sensi- 
tive to changes in modulus . Two more runs in the range of modulus 
indicated by the dotted lines would provide the information needed to 
construct the viscoelastic solution that will now be developed in detail. 




Crack Length (in) 

Fig. 23 

X vs . Crack Length 
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Consider first a stationary crack in this geometry. We have shown 
previously that the viscoelastic stress intensity factor for this case 
may be obtained from the correspondence principle (Eqn. 33). 

K l V s (t) = L" 1 [X*(E*,c) Og*(E*,P*) JT ] (65) 

V/ S 

a (t) can be found from the elastic solution by the method outlined in 
Appendix I. 

The viscoelastic correction factor, \ , can be found from Fig. 23 
by using an equivalent elastic modulus, which can be found for the viscoelastic 
problem. To illustrate, the equivalent modulus for plane strain is defined 
as follows : 



e 0 = IT^ C( 1_u ) ( 66a ) 

eq 

e 0 = C ( 1- u) Ug - uc/] (66b) 

eq 



In other words, the instantaneous viscoelastic stress and strain can be 
used to find the value of Young's modulus which would produce those 
stresses and strains if the body were elastic. Of course the equivalent 
modulus varies with time. 

The stresses and strains generally vary also from point to point 
within the motor; however for a cylinder the sum of normal stresses and 
strains is a constant throughout the motor. We will therefore choose to 
define the equivalent modulus in terms of this sum. Adding Eqns . (66a) 
and (66b) , 



v 

e + e . 
r e 



v _ (l+u) (l-2u) p 

° E L 

eq 



v v, 

<*- + CT J 



(67) 



42 



The viscoelastic strains can he found in terms of viscoelastic stresses 
through the viscoelastic constitutive law. 



®r + e e = ( 1+ u)(l-2u) J* D crp (t- T ) ^ (a^+apdT 



(68) 



0 - 



Eliminating (e r +e Q ) from Eqns. ( 67 ) and (68), the equivalent modulus 
becomes , 



E (t) -r- 
eq' ' t 



v v 

cr + 
r 9 



I D crp (t - T) 5F (ct ^ + CT 0 )dT 

0 - 



( 69 ) 



With this parameter determined, Fig. 23 can be used to find \^(t) = \[E e ^(t),c] 
An inversion of Eqn. ( 65 ) for a fixed crack length yields 



K 1>) = /c” J ^ v ( t_T ) <^(T)dT 

0 - 



(70) 



where 

\ v (t-r) = 0(E,c)] E = E eq (t-r) 
Referring back to Eqn. ( 69 ), 



E 



eq 



(t-T) 



aj(t-r) + aj(t-r) 

t-T 

£ D crp (t_T ' ?) 4 [a r (?) + CT I (5)]d§ 



(71) 



With the geometry and material correction factor included in the 
analysis so that a continuous updating of these properties can be effected, 



43 



the correlation between static and dynamic viscoelastic stress intensity 
factor would follow from Eqn. ( 38 ). 

t 

K^(t) = v^ItT J* [\(E,c)] c=c(t) ct q^( t )dr (72) 

0 E=E (t-T) 

Schematically the solution may he viewed to proceed as shown in Fig. 
24. The static viscoelastic solution is given for a variety of crack 
lengths. The viscoelastic stress intensity factor for a moving crack is 
assembled from the incremental pieces of the static solution. 




Fig. 2k 

Schmetic of Solution 
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Procedure 



To illustrate how the solution should proceed, a step by step set of 
instructions will be given. 

1. An appropriate time increment is selected and ct V and cf. are computed 

for p = p(t^) and E = from Eqn. (l-19) for any given 

pressure-time trace. (r is taken to be a + c^.) A time incre- 
ment is selected by guess using estimates of load and resulting 
crack velocities as guidelines. The time increment should allow 
for moderate extension of the crack but not be so large that the 
accuracy of the solution is compromised. 

2. Form a + o and compute the equivalent modulus from Eqn. ( 69 ) 

r 0 

for t = t^ and t = t^. In summation form it becomes, 



E (t ) 
eq v rr 



o V (t ) + oY(t ) 
r v n' 0 V n' 



n 

E 

i=0 



D (t -t.)[(a V (t.)+aJ(t.))-a V (t. ) - a J(t, -l)] 

crp n i /LV r v 1 0 1 r l-l 0 1 J 



3. Knowing E g ^'s for the first time increment, X^'s can be found for 
the corresponding times from Fig. 23. 

4. Knowing X from step 3 and aj from step 1, compute K 1 V (t 1 ) from 

"V 0 IS X 

Eqn. (70) for the first time increment. 



n 



=>n-l E h/W °e (t i )At i 
1=0 



5 . Using the fracture characterization curve for the grain propellant , 
(for example, see Fig. 17) find the crack velocity associated with 
the first time increment, c . Knowing the velocity and the length 
of the time increment, a new crack length, c^, can be computed to 
be used for the next iteration. 
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6. Progress to t = t g and 0 = 0 ^ and find cy^(tg) and , using 

p = p(tg), E = Ej^,^(tg) and r = a + for the second time 
increment and values as before for t^. 

7. Compute E g ^ for t = t^, t = t^, t = t Q . 

8. Find corresponding X^'s. 

9. Compute ^^(tg) and find the new crack velocity and crack length. 

10. The process is repeated until the entire crack history is 
recorded as shown in Fig. 24. 

11. From this information a crack history plot can be constructed as 
shown in Fig. 25. 



I 



■Quasi-viscoelastic Solution 



■ Viscoelastic Solution 



I 1 1 I 1 1 

tj. ^4 t 

Fig. 25 Viscoelastic Crack History 

This approximate solution is obtainable with realistic computer times 
and associated expenditures. At the same time it should yield good results. 
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SUMMARY AND CONCLUSIONS 



The basic definitions for viscoelastic stress intensity factors are 
set forth with associated computations for determining them in laboratory 
specimens. Laboratory tests have been conducted to characterize propellants 
fracturewise , and the most suitable test found was the long strip subjected 
to a constant force loading. The influence of dewetting on the fracture 
characterization appears to be significant, but the matter has not yet been 
fully resolved. 

An approximate method for performing viscoelastic fracture analyses 
has been presented which incorporates the fore -mentioned fracture character- 
ization. It requires minimum computer usage and should yield good results. 
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APPENDIX I 



An approximate method of viscoelastic analysis can be developed 
from linear superposition of elastic solutions. Such quasiviscoelastic 
solutions can give good solutions if properly used. The guidelines for 
use are formulated in this appendix by developing the quasiviscoelastic 
equations. 

The elastic normal stress expressions for an encased circular 
cylinder are of the form 

°r,e e - V + B l p t 1 - 1 * 

where p' is the interface pressure, and p is the bore pressure, and 
B^ are functions of motor geometry and the radial coordinate; however 
they are not time dependent. 

For , 




For cr 



0 



A l = 



-b 



b 2 - 



[1 



2 

a 

r 



B 1 - 2 



n if 
2 1 + 2 ] 
b - a r 



From the elastic analysis. 



P 



f 



P 

A 2 + b 2 e 



( 1 - 2 ) 



where A 2 and Bg are composed of geometrical constants, and E is Young's 
Modulus of the propellant. The normal stress expression can then be 
written as 



a 

r. 



_ P 
9e A + BE 



+ B-jP 



(1-3) 
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By the correspondence principle, the transformed stress for a unit 



step bore pressure is 



v* 



a = — iLB . — + JL 

r,0I A + BE* s 



(1-4) 

where the asterisk indicates the Laplace Transform and p* = l/s. The 

subscript I on stress indicates the indicial stress response; i.e. the 

stress resulting from a unit step input loading. 

For convenience in manipulation the first term in equation (1-2) 

* 

will be designated o^. . 



* 

CT 



111 



(1-5) 



I A + BE* 

Clearing Equation (l-5) of fractions and inverting to the time 
plane , 

t do 

n ( t) + B V o) E KEL (t) + B 4+ E m (t - T) "ai dT - 1 < x - 6) 

If the integral is broken up into time increments for evaluation by sum- 
mation, the expression (1-5) becomes 

n 



Aa l( t) + B CT ^( o) ERE1 (t) + B £ Erel (t - tl ) Aa n „ 1 (I- 



7) 



When t = t = 0, Equation (1-7) yields 



°i (o) ’ 



( 1 - 8 ) 



This is the same as Equation ( 1-3) with E replaced by E^^t). 
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Proceeding on in time to t = t^, Equation (l-7) now becomes, 

A VV + B V°> WV + B E BEL (o)to I (t l ) - 1 



where 



Aa( tl ) = - ^(o) 

If BO' T (t.) E n _ T (t 1 ) is added and subtracted to Equation (1-9), it can 

1 Kill Li _L 



be written in the form 



where 



[A + B E EEL (t l )] W - ^ = 1 f 1 ’ 10 ) 



E REl/M WV ■ E REL^ o) 



If A E^g L and Ao^ are chosen to be small compared to E^^ and in 
the numerical computing process, 



4 ’WV " A + B WV 



(I- 11) 



Once again this of the form of Equation ( I— 3 ) with E replaced by E OTT (t 1 ). 

Kill Li _L 

This procedure can be repeated for t = t^; however more manipulating is 
required. The result is 

f A * B E REL (t 2» W - B[4 E REL ( V 4 W + 4 E REL (t 2> 4 



+ AE HEL (t 2 )Aa I (t 2 )]=1 



( 1 - 12 ) 



By the same criterion as before, 



W 4 






(1-13) 
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It is assumed that this process can be repeated for an arbitrary 
number of time increments. As a result, a general relationship can be 
deduced 



CT ! (t i ) A + B E_ T (t.) 



( 1 - 14 ) 



KEL i J 

From equation (l-4) the viscoelastic stress components become upon 
inversion 



CT r,0I (t i ) ’ A + B E^ T (t.) + B 1 



KEL 



( 1 - 15 ) 



Proceeding back one step further, the viscoelastic stress field resulting 
from a unit step in bore pressure is found from equation (l-3). 



cr fl _(t.) 4 [ -1. 2*1 ] 

r,9l K P 



E = WV 



(1-16) 



Since 0I is the indicial response, the viscoelastic stress for 
an arbitrary pressure history can be obtained via the linear superposition 
integral. 




(t) 




CT r,ei^- T ) 



dp( T ) 

dT 



dT 



( 1 - 17 ) 



In terms suited to numerical computation. 



v 

a 

r 






e (t) = p(o)^ 0I (t) + 



n 

E 

i=l 



o a (t-t.)Ap. 
r,9l v i' 



(1-18) 
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Relating this back to the elastic stress expression through equation 



(1-16), 



a rj e (t) = p(o) 



°r , 9e n 

P E = E EEL (t) + 



n 

E 

i=l 



•,9e 



Pi 



(1-19) 



The viscoelastic stress is obtained by subdividing the time into increments, 
determining the associated increments in applied pressure, Ap^ , and then 
appropriately summing. 
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